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We derive the sub-subleading two-pion exchange contributions to the three-nucleon force which 
appear at next-to-next-to-next-to-next-to-leading order in chiral effective field theory. In order to 
determine the low-energy constants, a complete analysis of pion-nucleon scattering at the subleading- 
loop order in the heavy-baryon expansion is carried out utilizing the power counting scheme employed 
Cn . in the derivation of the nuclear forces. We discuss the convergence of the chiral expansion for this 

^"^ ' particular three-nucleon force topology and give the values of the low-energy constants which provide 

^^ ' the most realistic description of the three-nucleon force when the chiral expansion is truncated at 

Cn , next-to-next-to-leading order. 

C^ ■ PACS numbers: 13.75.Cs,21.30.-x 



I. INTRODUCTION 

I— { ' Three-nucleon forces (3NFs) are well known to play an important role in nuclear physics. In spite of many decades 
r^ ' of effort, the detailed structure of the 3NF is not properly reproduced by modern phenomenological 3NF models, see 
Ref. [l[ for a comprehensive review. This provides a strong motivation to explore the structure of the 3NF within chiral 
effective field theory (EFT) which is nowadays a standard tool to analyze low-energy nuclear dynamics in harmony 
with the symmetries of QCD, see Refs. [2,-^4^] for recent review articles. In the past two decades, nuclear forces have 
r~^ , already been extensively studied in this framework. For two nucleons, it turned out to be necessary and sufficient to 
\^ ' go to next-to-next-to-next-to-leading order (N'^LO) in order to accurately describe the static deuteron properties as 
^^ , well as the two-nuclcon phase shifts and mixing angles up to laboratory energies of i^iab ^ 200 MeV [g, Q • On the 
^D ■ other hand, three- and more-nucleon systems are so far only analyzed up to next-to-next-to-leading order (N^LO) 
^f^ ' in the chiral expansion [3, |8[ . At this order one gets the first non- vanishing contributions to the 3NF which emerge 
f^ from the two-pion-exchange, one-pion-exchange-contact and purely contact graphs (a), (d) and (f) in Fig. [T] with 
CN [ the corresponding amplitudes being given by the lowest-order pion-nucleon vertices. Generally, one observes a good 
description of nucleon-deuteron elastic and breakup scattering observables at very low energies which improves when 
going from next-to-leading order (NLO) to N^LO. On the other hand, the well-known puzzles in the three-nucleon 
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continuum such as e.g. the ^j^-puzzle 



HI 



and the large discrepancy for the breakup cross section in the so-called 



space-star and related configurations [l^ll3| still persist at N^LO. Notice, however, that the recent calculation by the 

j^ [ Pisa group [ij] demonstrates that the ^y-puzzle in the 4N system is significantly reduced by the chiral 3NF at N^LO. 

" " The chiral EFT predictions for the three-nucleon scattering observables at N^LO at intermediate and higher energies 

are, in general, in agreement with the data but show a rapidly increasing theoretical uncertainty. It is, therefore, 

necessary to go to higher orders in the chiral expansion for three- and more-nucleon systems. 

The N'^LO contributions to the 3NF emerge from the leading relativistic corrections and pion-loop diagrams in all 
six topologies shown in Fig. [TJ It is important to stress that the N^LO contributions do not involve any unknown 
low-energ y c onstants (LECs). The corresponding parameter- free expressions can be found in Refs. [H, [l^, see 
also Ref. [I4I . Another interesting feature of the N^LO 3NF corrections is their rather rich isospin-spin-momentum 
structure emerging, especially, from the ring topology (c) in Fig.[T] This is in contrast with the quite restricted operator 
structure of the two-pion exchange 3NF topology (a) whose effects in the three-nucleon continuum have already been 
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FIG. 1: Various topologies contributing to the 3NF up to and including N^LO: two-pion {2-k) exchange (a), two-pion-one-pion 
(27r-l7r) exchange (b), ring (c), one-pion-exchange-contact (d), two-pion-exchange-contact (e) and purely contact (f) diagrams. 
Solid and dashed lines represent nucleons and pions, respectively. Shaded blobs represent the corresponding amplitudes. 



extensively explored. It is, therefore, very interesting to study the impact of the novel structures in the 3NF on 
nucleon-deuteron scattering and the properties of light nuclei, especially in connection with the already mentioned 
unsolved puzzles. On the other hand, one may ask whether the resulting (leading) contributions to the structure 
functions accompanying the novel operator structures in the 3NF already allow for their decent description. Stated 
differently, the question is whether the lowest-nonvanishing-order contributions from the 27r-l7r and ring-topologies are 
already converged or, at least, provide a reasonable approximation to the converged result. There is a strong reason to 
believe that this is not going to be the case since the contributions due to intermediate A(1232) excitations are not yet 
taken into account for these topologies at N'^LO. In the standard chiral EFT formulation based on pions and nucleons 
as the only explicit degrees of freedom, effects of the A (and heavier resonances as well as heavy mesons) are hidden 
in the (renormalized) values of certain LECs starting from the subleading effective Lagrangian. The major part of 
the A contributions to the nuclear forces is taken into account in the A-less theory through resonance saturation of 
the LECs C34 accompanying the subleading nnNN vertices [18l - (22 | (see, however, the last two references for some 
examples of the A-contributions that go beyond the saturation of 03^4). These LECs turn out to be numerically large 
and are known to be driven by the A isobar [20,l23|. As a consequence, one observes a rather unnatural convergence 
pattern in the chiral expansion of the two-pion exchange nucleon-nucleon potential V^Z with by far the strongest 
contribution resulting from the formally subleading triangle diagram proportional to C3 |2j|. The (formally) leading 
contribution to V^^ does not provide a good approximation to the potential so that one needs to go to (at least) the 
next-higher order in the chiral expansion and/or to include the A isobar as an explicit degree of freedom [20| . The 
situation with the 2tt-1tt and ring topologies in the 3NF is similar. Based on the experience with the two-nucleon 
potential, one expects significant contributions due to intermediate A excitations, see also the discussion in Ref. ,25.] . 
For the ring topology, this expectation is confirmed by the phenomenological study of Ref. [20| . In order to include 
effects of the A-isobar one needs 



• either to go to (at least) next-to-next-to-next-to-next-to-leading order (N^'LO) in the standard A-less EFT 
approach 

• or to include the A-isobar as an explicit degree of freedom. 



It should be understood that both strategies outlined above are, to some extent, complementary to each other. In 
particular, N'^LO contributions in the A-less theory only take into account effects due to single A-excitation but not 
due to the double and triple A-excitations (whose inclusion in the A-less approach would require the calculation at 
even higher orders). These effects are taken into account already at N'^LO in the A-fuU approach. On the other 
hand, there are also contributions not related to A-excitations which are included/absent in the A-less approach at 
N^LO/A-fuU theory at N^LO. It remains to be seen which strategy will turn out to be most efficient. The present 
paper represents the first step along this line. We analyze here the longest-range contribution to the 3NF in the 
standard, A-less approach at N'^LO in the chiral expansion. This topology is particularly challenging due to (i) the 
need to carry out a non-trivial renormalization program as explained in section IIIII and (ii) the need to re-consider 
pion-nuclcon scattering in order to determine the relevant LECs. Our paper is organized as follows. In section |lll we 
specify all terms in the effective Lagrangian that are needed in the calculation. The general structure of the two-pion 
exchange 3NF is discussed in section [1111 Here, we also briefly summarize the already available results at N^LO and 
N'^LO and give explicit expressions for the N^LO contributions. In section HVl we analyze pion-nucleon scattering at 
order Q^ in the chiral expansion with Q referring to the soft scale of the order of the pion mass and use the available 
partial wave analyses to determine the relevant LECs. In section |Vl the numerical results for the two-pion exchange 



3NF are presented and the convergence of the chiral expansion is discussed. Finally, the findings of our work are 
briefly summarized in section IVII The Appendices contain explicit formulae for the algebraic structure of the nuclear 
Hamiltonian at order N'*LO and the chiral expansion of the ttN invariant amplitudes. 



II. EFFECTIVE LAGRANGIAN 



To derive the longest-range contributions to the 3NF at N'^LO we need the effective pion-nucleon Lagrangian up to 
the order Q^ . The explicit expressions in the heavy-baryon formulation can be found in Ref . [27|, 123 • I^O'^ ^^^ sake of 
completeness, we list here all terms relevant for our calculation with the corresponding building blocks being expanded 
in powers of the pion fields: 
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where tt and Ny refer to pion and nucleon fields, t denote the isospin Pauli matrices, v is the nucleon four-velocity 
and Sf^ refers to the covariant spin operator of the nucleon. 
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Further, F and qa are the pion decay and the nucleon axial vector constants in the chiral limit, M is the pion mass 
to leading order in quark masses and di (e^) are further low-energy constants (LECs) from the order-Q'^ (order-Q'') 
pion-nucleon Lagrangian. The superscript i of C^j^ refers to the number of derivatives or insertions of the pion mass. 
Notice that in the power counting scheme employed here, the nucleon mass is treated as a heavier scale as compared 
to the chiral-symmetry-breaking scale A^^, see the discussion at the end of the section UlI Bl As a consequence, the 

l/m2-terms in £^^ and the l/m-, Xjm}- and l/m'^-terms in O-^]^ generate contributions to the 3NF beyond N^LO. 



We, therefore, refrain from listing these terms in the effective Lagrangian. The constant a represents the freedom in 
parametrizing the pion field. Clearly, physical observables should not depend on a particular choice of a. Therefore, 
keeping a. unspecified and verifying a-independence of the resulting nuclear forces allows for a non-trivial check of 
the calculation. Notice that isospin-breaking corrections to the 3NF start to contribute at N^LO. The complete 
expressions for isospin-breaking terms in the 3NF up to and including N^'LO can be found in Ref. [29[, see also 
[30, |31|. We, therefore, refrain from including isospin-breaking terms in the present work and employ exact isospin 
symmetry. For more details on the effective pion-nucleon Lagrangian the reader is referred to Refs. [28l l32l l33l]. 



III. THE TWO-PION EXCHANGE 3NF 

The 27r-exchange topology (a) generates the longest-range contribution to the 3NF. In the isospin and static limits, 
its general structure in momentum space has the following form (modulo terms of a shorter range such as e.g. the 
ones corresponding to the (b)-topology): 

^'" " {q\ + M^] S -F M2] (^' ■ ^^ -^^'^'^ + ^1 X T3 • X2 qi X 93 • (72 6(92)) , (3.3) 

where Oi denote the Pauli spin matrices for the nucleon i and qi — Pi' — Pi, with pi' and pi being the final and 
initial momenta of the nucleon i. Here and in what follows, we use the notation: qi = \qi\. The quantities A{q2) and 
B{q2) in Eq. p.3p are scalar functions of the momentum transfer (72 of the second nucleon whose explicite form is 
derived within the chiral expansion. Unless stated otherwise, the expressions for the 3NF results are always given for 
a particular choice of the nucleon labels. The complete result can then be found by taking into account all possible 
permutations of the nucleons 

^3N ' = ^3N + 5 permutations . (3.4) 

It is important to emphasize that we are using in Eq. p.3p a slightly different notation as compared to e.g. Ref. |15| . 
Specifically, in order to avoid the issue of non-uniqueness of a decomposition into the 2-k, 2tt-1tt and shorter-range 
contributions, we get rid of all terms which are proportional to qf and gf. More precisely, using the identity (?i 3 = 
{qf 3 -I- M^) — M^ and canceling the inverse pion propagator with the corresponding one in Eq. p.3p . each qf 3 in the 
numerator gets replaced by —M^ modulo some additional contributions to the two-pion-one-pion exchange, one-pion- 
exchange-contact and the purely short-range contact interactions Kont , see graphs (b) , (d) and (f ) in Fig. [1] This way 
we ensure that the resulting functions A and B depend solely on (72 rather than on qi, (72 and q^. The purely short- 
range, induced terms T4ont only shift the low-energy constants accompanying the contact 3NFs whose values anyway 
need to be adjusted to the data. There is no need to keep these contributions explicitly. The induced contributions 
to the other two topologies do, however, need to be taken into account, see Ref. [3J| for a related discussion. 



A. N^LO and N^LO contributions 

We now briefly consider the first two terms in the chiral expansion of the functions A{q2) and B{q2). The leading 
contributions arise at N^LO which corresponds to the order Q^ relative to the leading contribution to the nuclear 
Hamiltonian and have the form [3, Q 

A^'Hq2) = ^ ((2C3 - ic,)AC + c^ql) , B^'Hq2) = ^ , (3.5) 

where qa, F^^ and M^ denote to the physical values of the nucleon axial vector coupling, pion decay constant and 
pion mass, respectively, and the superscripts of A and B refer to the powers of the soft scale Q. The first corrections 
at N^LOread ^\l^: 

l3^'Hq2) - - 2^^ [^(92) (4M^ + 9I) + (23^ + mu] , (3.6) 
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FIG. 2: Two-pion exchange 3N diagrams at N*LO. Solid dots, filled circles, rectangles and diamonds denote vertices from 
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C^ ' or Z"-" 
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, respectively. Open rectangles (diamonds) refer to l/77i-vertices from C^l^ 
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Diagrams which result from the interchange of the nucleon lines and/or application of the time reversal operation are not 
shown. For remaining notation see Fig. [l] 



where the loop function A{q) is defined as: 



A(q) — —- arctan -— -— 



(3.7) 



Notice that the leading-loop contributions to the 27r-exchange topology do not contain logarithmic ultraviolet diver- 
gences and, as explained in Ref. US'], turn out to be independent from the LECs di entering C^^lf. At both N^LO 
and N'^LO, all LECs in the effective Lagrangian entering the expressions for the 3NF - including gA and Ft^ - can be 
simply replaced by their physical values. 

We further emphasize that, as already mentioned above, the expressions in Eq. p.6p differ from the ones in Eq. (2.9) 
of Ref. [151 by terms of a shorter range as compared to the two-pion exchange contributions. The advantage of using 
the new notation is that the results for A and B are now a-independent. This was not the case for terms in Eq. (2.9) 
of Ref. [l^ where the results are given for a specific choice a = 0. 

Last but not least, we emphasize that relativistic corrections to V27r have a richer structure than the one given in 
Eq. p.3p . The explicit form of the 1/m-corrcctions to V27r at N^LO can be found in Ref. [ly|, see also [3a] for an early 
work. 



B. N LO contributions 



We now turn to the sub-subleading contributions to the 27r-exchange 3NF at order Q^ (N^LO). These are depicted 
in Fig. [5] and emerge from: 



• one-loop diagrams (1)-(15) constructed from the leading-order vertices from C^jf and a single insertion of a 
subleading vertex c>c Ci, 

• tree diagram (16) involving leading-order vertices from /Cj|.jv ^^"^ ^ single insertion of a vertex from >C^^ pro- 
portional to LECs Bi, 



• tree diagram (17) involving leading-order vertices from C)^^, one £^^-vertex oc q and one £^^-vertex oc di, 

• relativistic l/rn-corrections resulting from diagrams (18)-(21). 

Notice that since we work with renormalized pion field operators, we do not show explicitly in Fig.[2]diagrams involving 
pion self energy. We further emphasize that some of the diagrams shown in Fig. [2l such as e.g. graphs (10) and (11), 
yield vanishing contributions to the 3NF. 

It is important to keep in mind that, in order to derive the genuine 3NF contributions, we need to separate the 
irreducible parts in the corresponding amplitudes to avoid double counting when iterating the potentials in the 
scattering equation. This is achieved employin g th e method of unitary transformation, see Refs. [36i, l37J for a 
comprehensive description of the method, Refs. [H, [Ifl, |33, |33 for recent higher-order calculations of three- and 
four-nuclear forces and Refs. [40, |41| for an extension to electromagnetic processes. More precisely, we only use 
the method of unitary transformation to evaluate the contributions of graphs (1), (3), (5), (6)-(8), (20) and (21) in 
Fig. [T] which involve reducible topologies. The remaining contributions are obtained by calculating the corresponding 
Feynmann diagrams. 

Before showing our results for the N^LO contributions to the functions A and B^ we briefly remind the reader how the 
calculations within the method of unitary transformation are organized. We begin with the effective Lagrangian spec- 
ified in section [n] and switch to the Hamiltonian using the canonical formalism. We prefer to work with renormalized 

— 1/2 

pion field tt*" = Z^ ' tt and mass M^ and, therefore, do not need to deel with pion self-energy contributions. In the 
second step, the pion degrees of freedom are projected out by employing the appropriate unitary transformation in 
the Fock space. The explicit form of the unitary operator needed to compute nuclear forces up to N^LO can be found 
in Ref. [Ig, [S^ • It is important to emphasize that getting rid of ultraviolet divergences appearing in the pion loop 
contributions to the nuclear Hamiltonian and current operators by expressing the bare LECs in terms of renormalized 
ones is not guaranteed a priori (given that the nuclear forces and currents are not observable quantities). This issue is 
discussed in detail in Ref. 139| , see also Ref. [4l| . To ensure renormalizability, we exploit the unitary ambiguity of the 
resulting nuclear potentials, i.e. the freedom in choosing the basis states in the Hilbert space. More precisely, we apply 
additional unitary transformations in the purely nucleonic subspace of the Fock space whose "rotation angles" are 
chosen in such a way, that all ultraviolet divergences are absorbed into redefinition of the LECs. As explicitly demon- 
strated in Refs. [3^, |4l|, it is possible to carry out this renormalization program for the leading-loop contributions to 
the 3NF and the electromagnetic two-nucleon current operators. An extension to subleading-loop 3NF contributions 
at N'*LO is straightforward. In Appendix |X] we show the resulting formal algebraic structure of the N'*LO contri- 
butions to the nuclear force ex g\ci and the retardation corrections ex g\ci/m. To evaluate the corresponding 3NF 
contributions, we simply compute the matrix elements of the connected time-ordered-like 3N diagrams emerging from 
these operators and employ dimensional regularization to deal with ultraviolet divergencies. The unitary ambiguity is 
parametrized via the "rotation angles" ag.io.ii which have to be chosen in such a way that the 3NF matrix elements 
become finite when expressed in terms of renormalized LECs. At the order we are working, the relations between the 
LECs (jA, F and M and the corresponding renormalized constants are given by 
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where Z34 are the LECs entering the subleading pion Lagrangian Ctt |42,] and the (divergent) quantity A^r is defined, 
following the notation of Ref. [4J] , via 
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We also need the expression for the pion Z-factor Z^r, which has the form 

2A/2 



(3.9) 



^- = 1 + ^ ((10" - 1)^- - ^4) + 0{Mt) . (3.10) 



Notice that the Z-factor is not observable and shows an exphcit dependence on a. For a = 0, our result coincide 
with the one given in Ref. 4J]. Further, the LECs li, di and e^ can be decomposed into the divergent parts and finite 



pieces. Utilizing the notation of Ref. [43]^ this decomposition has the form: 

Ci = e[(^) + -^A = e, + -^A^, (3.11) 

where /i denotes the renormalization scale and the divergent, /i-dependent quantity A is related to the /^-independent 
one Ajr through 

A^A.-^log(-). (3.12) 

The coefficients % and Si are well known in the framework of dimensional regularization adopted in the present work 
[Ullllillill. The relevant coefficients read: 

Sis = 0, Sie = -gA + gA, 73 = -^, 74 = 2. (3.13) 

To the best of our knowledge, the ei-coefficients have not yet been worked out. Fortunately, we only need a few linear 
combinations of e^ which are exactly the same as ones appearing in ttN scattering at order Q'*. Cancellation of the 
ultraviolet divergences in the irN scattering amplitude implies the following relations for the divergent parts of the 
LECs e, HI: 

1 1 

5 3 

-2ei9 + £22 + £36 = -2ci + — C2 - -C3 , 

2 2 

-2e2i + £37 = -o'^4-3g^C4, 
1 3 

£22 - 4638 = -3ci + -C2 + -C3 . (3.14) 

Expressing the N'^LO contribution to the 3NF in terms of gA, Mt^, Fj^ and the renormalized LECs k, di, e^ using 
the above relations leads to finite matrix elements provided the "rotation angles" ag_io,ii of the additional unitary 
transformations are chosen as: 

aio = --(1 -2a9), an = -(l-2a9). (3.15) 

Notice that while the parameter ag is unfixed, the resulting 3NF turns out to be ag-independent. This leads to an 
unambiguous result for the 3NF at this order. We further emphasize that the obtained relations for aio,ii also constrain 
the form of the remaining N^'LO 3NF contributions, the four-nucleon force and the subleading-loop expressions for 
the exchange current operators. 

The final, renormalized N^'LO contributions to the functions A and B in Eq. p.3p have the form: 

A^'H<12) = ^^^|^[Af^(7l(i^'(2304^V(4ei4 + 2ei9-e22-e36)-23047r2Ji8C3) 

+ .gA(144ci - 53c2 - 9OC3)) + Af^ {F^ {A608n^dis{2ci - C3) + A608n^gA{2eu + 2ei9 - 636 - 4e38)) 



^ To simplify tlie notation, we use here diniensionful e^ in contrast to that reference. 
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where the loop function L(q) is defined according to 
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Interestingly, we observe that there are no 1/m-contributions to the two-pion exchange 3NF at this order. In particular, 
no N^LO contributions emerge from diagrams (18) and (19) in Fig.[5]since all leading (subleading) 1/m-corrections to 
the ttNN (jTirNN) vertex involve at least one time derivative. When evaluating the corresponding Feynman diagrams, 
these time derivatives generate insertions of the nucleon kinetic energy which are further suppressed by the factor 
Q/m. In addition, diagrms (20) and (21) are found not to generate any irreducible pieces. 

The expressions for the 27r-exchange 3NF up to N'^LO discussed above depend on a number of low-energy constants. 
Here and in what follows, we use the values^ 



5A = 1.267, i^^=92.4MeV, 



M^ = 138.03 MeV. 



(3.18) 



The LECs Ci, di and e^ can be most naturally determined from pion- nucleon scattering (at least) at the sublcading- 
loop order (i.e. Q^). The heavy-baryon analyses of pion- nucleon scattering at orders Q^ and Q'^ can be found in 
Refs. [32, Sa |46[ , see also Refs. [4J, |47|, |48J for the calculations within the manifestly covariant framework, Ref. [4^ 
for a related calculation which extends chiral EFT to higher energies by employing constraints set by causality and 
unitarity and Ref. [50| for a recent review on baryon chiral perturbation theory. Unfortunately, we cannot use the 
values of the LECs obtained in these studies since we use a different counting scheme for the nucleon mass in the 
few-nucleon sector, namely Q/m ~ Q'^ / h.^ jSlj] rather then m ~ A^ as used in the single-nucleon sector, see [S^l for 
an extended discussion. In the next section, we re-analyze pion-nucleon scattering at order Q in the heavy-baryon 
approach utilizing our counting scheme for the nucleon mass and determine all relevant LECs from a fit to the available 
partial wave analyses. 



IV. DETERMINATION OF THE LECS FROM -kN SCATTERING AT ORDER Q* 



In the center-of-mass system (cms), the amplitude for the reaction ■n'^iqi) -\- N{pi) — ^ 71^(^2) + N{p2) with pi^2 and 
91,2 being the corresponding four-momenta and a, b referring to the pion isospin quantum numbers, takes the form: 



rpba 



E 



2ni 



Tba 



g^{u},t) + ia ■ q2 X qi h'^{u},t) 



ie'""=T'^ 



g {uj,t) +i<T ■ q2 X qih {uj,t) 



(4.19) 



E 

^2 



(q^ + m2)^/2 ^jjg nucleon energy and qi ^ 



92' 



Here, uj = qi = q2 is the pion cms energy, Ei = E2 

{{s — M^ — m^)^ — Am? M"^) / {As) . Further, t = [qi — 92)'^ is the invariant momentum transfer squared while s denotes 
the total cms energy squared. The quantities g {uj,t) {h {uj,t)) refer to the isoscalar and isovector non-spin-fiip 
(spin-fiip) amplitudes and can be calculated in chiral perturbation theory. In Appendix [B] we show the contributions 
to the amplitudes up to and including the order Q^ using the same counting scheme for the nucleon mass as in the 



^ Since we employ exact isospin limit in this work, we do not distinguish between the charge and neutral pion masses and use M^ 

2M^+/3 + M^o/3. 



derivation of the nuclear forces. It is important to emphasize that terms in Eq. (IB.4[) proportional to the LECs eig, 
e20j e2ij e22j 635, 630, 637, Cag and I3 can be generated from the Q^-terms by shifting the LECs q as follows j45| : 



Ci ^ Ci - 2Af^ £22 - 4638 






C2 -^ C2 + 8Af2(e2o + e35), 



C3 ^ C3 + 4M^ ( 2ei9 - 622 - 636 + 2-^2 , , 

C4 -^ C4 + 4Af2(2e2i-e37). (4.20) 

This implies that one cannot extract these combinations of Ei and Ci separately from the ttN scattering data. We, 
therefore, follow Ref. [4^ and absorb these linear combinations into redefinition of the c^'s by setting 

l^Ci _ _ _ _ _ ^361 

622 - 4638 ~ -^ = 0, 620 + 635 = 0, 2ei9 " 622 " 636 + 2—^ = 0, 2621 - 637 = , (4-21) 

^ IT -^ 7T 

without loss of generality. Choosing another prescription would result in corrections generated by the loops involving 
Ci that are beyond the order we are working. Notice further that with the above convention, there is no dependence 
anymore on the LEC ^3. The LEC rfis can be fixed by means of the Goldberber-Treiman discrepancy 

Ftt V 9A ) 

where for g-irNN we adopt the value from Ref. |53i] : 5^ArAr/(47r) ~ 13.54 which also agrees with the recent determination 
in Ref. [5J| based on the Goldberger-Miyazawa-Oehme sum rule and utilizing the most accurate available data on the 
pion-nucleon scattering lengths. Again, at the order we are working, we are free to set dis — provided we use the 
effective value for gA, 

^ F^g^NN ^ ^285, (4.23) 

m 

in all expressions. We adopt this convention in the following analysis. Thus, we are finally left with 13 independent 

(linear combinations of the) low energy constants to be fixed from a fit to the data, namely ci^2,3,4, di + (i2, (^3, ^5, 

dii — di5 and 614^15^16,17,18- 

The fit can be most conveniently performed in the partial wave basis using the available partial wave analyses. In 
order to estimate a possible uncertainty of the extracted parameters, we considered two different partial wave analyses 
in our fitting procedure, namely the one of Ref. [56| by the group at the George Washington University, to be referred 
as GW, and the Karlsruhe-Helsinki analysis of Ref. [53]), to be referred as KH. The energy range of the data fitted 
corresponds to the ttN laboratory momenta pLab < 150 MeV/c. At higher energies the convergence of the chiral 
expansion becomes doubtful. We follow the strategy which is similar to the one utilized in Ref. ^] and assign the 
same relative error to all empirical data points. 

The partial wave amplitudes /;_|_(s), where I refers to the orbital angular momentum and the subscript '±' to the 
total angular momentum (j = Z ± s), are given in terms of the invariant amplitudes via 



± , , E + m f+^ 



ft{^)-^^F^ d^ 



g^ Pi{z) + r h^ {Pi±i{z) - zPi{z)) 



where z ~ coii[9) is the angular variable (t = 2(f^(z — 1)). The amplitude in the isospin basis are related to /^^ as 
follows 

fi± = h± + 2//± , fi± = fi± - fi± ■ 

The phase shift for a partial wave amplitude with isospin / is obtained using the following unitarization prescription'^ 

5l^ {s) = arctan {\q\ 3? //± (s)) , (4.24) 



^ It should be understood that this unitarization prescription goes, strictly speaking, beyond the chiral power counting. The resulting 
model dependence is, however, very small due to the smallness of the phase shifts with the only exception of the P33 partial wave, see 
I55II for a related discussion. 
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FIG. 3: Results of the fit for tyN s, p and d-wave phase shifts using the GW partial wave analysis of Ref. j5q]. The solid curves 
correspond to the full Q results, the dashed curves to the order-Q results, and the dashed-dotted curves to the order-Q 
calculation. 

which reflects the absence of inelasticity below the two-pion production threshold. 

We performed a combined fit for all s-, p-, and d- waves since d- waves are the highest partial waves where the order-Q^ 
counter terms contribute. The results of the fits using the GW and KH partial wave analyses are visualized in Figs. [3] 
and m respectively. In these figures we show the full, order-Q'' results (solid curves) as well as the phase shifts 
calculated up to the order Q^ (dashed curves) and Q^ (dashed-dotted curves) using the same parameters (from the 
order-Q^ flt) in all curves. In the fitted region (from threshold up to pLab = 150 MeV/c), a good description of the 
data is achieved. As one would expect the convergence pattern when going from Q^ to Q^ is getting worse with 
increasing the pion momenta. Interestingly, the d-waves are rather well reproduced already at the order Q^ where 
there are no counter terms or other contributions depending on free parameters. Both the tree-level and finite loop 
contributions are important for those four partial waves. Our results for the phase shifts are similar and of a similar 
quality as the ones reported in Ref. [45[. 

We flnally turn to the discussion of the extracted parameters. The obtained values of the low energy constants are 
collected in Table HI As one can see from the table, the LECs Ci and di turn out to come out rather similar for the 
two partial wave analyses. The difference does not exceed 30% except for the LECs Ci and d^ which are, however, 
considerably smaller than the other Ci's and d^'s, respectively. The same conclusion about stability can be drawn for 
the LECs ei4 and 617. These are the only counter terms contributing to d- waves, which is why these two constants 
are strongly constrained by the threshold behavior of the d-wave phase shifts. In contrast, the other e^'s are very 
sensitive to the energy dependence of the s- and p-wave amplitudes and, therefore, vary strongly from one analysis to 
another. Notice, however, that all extracted constants are of a natural size except for the combination di4 — dis and 
ei5 which appear to be somewhat large. 



We stress that one cannot directly compare the LECs di and e^ from of our fits to the ones obtained in Refs. 
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FIG. 4: Results of the fit for tvN s, p and d-wave phase shifts using the KH partial wave analysis of Ref. p^. The solid curves 
correspond to the full Q results, the dashed curves to the order-Q results, and the dashed-dotted curves to the order-Q 
calculation. 

using heavy-baryon chiral perturbation theory at orders Q'^ and Q*, respectively, because of a different power counting 
schemes in the two approaches. On the other hand, it is comforting to see that the extracted values for the c^-, di- 
and even some of the Ci-coefficients are comparable to the ones found in Ref. [4^ in the fit with the LECs c^ being 
set to their order-Q^ values, see table 4 of that work. We also stress that the values for Ci.34 obtained from the fit 
to the KH partial wave analysis are in an excellent agreement with the ones determined at order Q^ by using chiral 
perturbation theory inside the Mandelstam triangle [58| . It is also worth mentioning that the values of 03^4 are in a 
good agreement with the ones determined from the new partial wave analysis of proton-proton and neutron-proton 
scattering data of Ref. [53 . 

It should be emphasized that one can obtain a considerably better description of the irN phase shifts at orders Q^ 
and Q^ by allowing for the LECs Ci and di to be tuned rather than keeping their values fixed at order Q^. In fact, 
the values of Ci are well known to change significantly when performing fits at orders Q^ and Q^. Using the KH 
partial wave analysis, employing the order-Q'^ expressions for the amplitudes and utilizing the same fitting procedure 





Cl 


C2 


c-i 


C4 


di +d2 


da 


ds 


di4 — di5 


ei4 


ei5 


ei6 


ei7 


ei8 


fit to GW, Ref. [56] 


-1.13 


3.69 


-5.51 


3.71 


5.57 


-5.35 


0.02 


-10.26 


1.75 


-5.80 


1.76 


-0.58 


0.96 


fit to KH, Ref. [57] 


-0.75 


3.49 


-4.77 


3.34 


6.21 


-6.83 


0.78 


-12.02 


1.52 


-10.41 


6.08 


-0.37 


3.26 



TABLE I: Low-energy constants obtained from a fit to the empirical s, p- and d-wave pion-nucleon phase shifts using partial 
wave analysis of Ref. [5y] and of Ref. [57|]. Values of the LECs are given in GeV^^, GeV~^ and GeV^'' for the d, di and en, 
respectively. 
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as before, we end up with the following values for the c^'s: 

cf" = -0.26 GeV"\ c^" = 2.02 GeV"\ c^" = -2.80 GeV-\ cf " = 2.01 GeV-\ (4.25) 

while the GW partial wave analysis yields: 

cf^ = -0.58 GeV-\ c^^ = 2.02 GcV-\ c^^ = -3.14 GeV-\ c^^ = 2.19 GeV'^ (4.26) 

Notice that 02,3,4 turn out to be somewhat smaller in magnitude than the ones extracted from the order-Q^ fit to the 
s- and p-wave TriV threshold coefficients ^J.'' We will come back to the issue of optimizing the description of the 
data at lower orders of the chiral expansion in the next section. Notice, however, that in spite of such a possibility, we 
believe that the results shown in Figs. |3] and |3] provide a more realistic picture of the convergence of the heavy-baryon 
chiral perturbation theory for pion-nucleon scattering. 



V. RESULTS FOR THE TWO-PION EXCHANGE 3NF 

With all relevant LECs being determined from pion-nucleon scattering, we are now in the position to analyze the 
convergence of the chiral expansion for the two-pion exchange 3NF. In Fig. [SJ we show the results for the functions 
A{q2) and B{q2) for small values of the momentum transfer (72, 92 < 300 MeV at various orders in the chiral expansion. 
More precisely, we plot ^(3), ^(3) + ^(4) ^nd ^f^) + ^(4) ^_ _4(5) ^^g ^g^ ^^g g(3)^ g(3) _^ ^(4) ^^^ ^(3) _^ g(4) _^ ^(5) 

using the values of the LECs q, di and e^ determined from the order-Q* fit to the KH and GW partial wave analyses 
as described in the previous section. We use here the same, fixed values for the LECs Ci (and di) listed in Table 
H] at all orders and adopt the same conventions regarding the LECs as in the case of pion-nucleon scattering, see 
Eqs. (|4.21l) and (I4.23p . Notice that A^^^ and B^^'' do not depend on the LECs 615,16,18 which are very sensitive to a 
particular choice of the partial wave analysis in pion-nucleon scattering, see Table HI The relevant LECs 614,17 are, on 
the contrary, rather stable as they are well determined in t:N d-waves. 

One observes a very good convergence for the function A with the subleading-order result (i.e. N^LO) being very 
close to the one at N'^LO. It is also comforting to see that both partial wave analyses lead to similar results for this 
quantity. The dependence on the input for pion-nucleon phase shifts for A is bigger than the shift from N^LO to 
N^'LO which can serve as a (conservative) estimation of the theoretical uncertainty at N'^LO. The convergence for the 
function B is somewhat slower with the shift from N^LO to N^'LO being of the order of ^ 30%. Also the difference 
between the two partial wave analyses of the order of ~ 20% is larger than for the function A. It should be understood 
that an accurate description of the low-energy pion-nucleon scattering data at different orders does not automatically 
guarantee a good convergence of the chiral expansion for A and B. In particular, these quantities do not depend on the 
LECs di (to the order considered) which contribute to nN phase shifts. Thus, the observed reasonable convergence 
for the 27r-exchange 3NF is a highly non-trivial test of the theoretical approach. 

Given that the (72 -dependence of A, B does not change significantly when going from N^LO to N'^LO, it is clear that 
the final result at N'^LO can be well approximated by the N^LO expressions with the appropriately tuned LECs 61,3,4. 
This feature is visualized in Fig. [HI where the dotted lines show .4^3) (g^) and S^^)^^^) from Eq. (|3.5p with 

cfH ^ _0.37 GeV-\ e™ = -2.71 GeV~\ 6™ = 1.41 GeV"^ (5.27) 

These values of o^'s are adjusted in such a way that yl(3)(0), B''^\Qi) and the curvature of the function A'^^\q2) at 
92 = coincide with the ones resufiing from y^^^) ^ _4(4) _|_ j^(b) ^ ^(3) _^ ^(4) ^_ ^(5) .^^^^ ^gj^^g ^j^g \MCs from the KH 

fits to the -kN data. Not surprisingly, the Oj's in Eq. (|5.27p are fairly close to the Q^-fit values in Eq. (|4.25p . We also 
show in Fig. inithe functions A^^\q2), B^^\q2) with the 6i's from Eq. (I4.25P (dashed-dotted lines). When using the 
GW partial wave analysis, the resulting values of the 6i's slightly change to: 

efw = -0.73 GeV^\ e^^ = -3.38 GeV~\ e^^ = 1.69 GeV-^ (5.28) 

The values in Eqs. ()5.27p and (|5.28p can be regarded as the recommended ones when the 3NF is taken into account 
at N^LO. One should, however, always keep in mind that the true theoretical uncertainty for the 27r-exchange 3NF 
at N^LO is sizable, see Fig. [5] 



* This indicates that the order-Q^ representation of the amplitudes does not provide the appropriate desciption of the data in the whole 
momentum range used in our fits. 
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FIG. 5: Chiral expansion of the functions A{q2) and B{q2) entering the two-pion exchange 3NF in Eq. (|3.3p up to N LO. Left 
(right) panel shows the resuhs obtained with the LECs determined from the order-Q* fit to the pion-nucleon partial wave 
analysis of Ref. [57| (Ref. [SQ]) and listed in Table|Tl Dashed, dashed-dotted and solid lines correspond to A''^' , A^^' + A^'^' and 

_4(3) ^ ^(4) ^ J^{5) jjj ^jjg yppgj. pjQ^g ^jjjjg ^(3)^ ^(3) ^ ^(4) ^^^ ^(3) ^ g(4) ^ ^(5) -^^ ^-^^ ^OWeV plots. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we have analyzed the longest-range contribution to the three-nucleon force at N*LO utilizing the heavy- 
baryon formulation of chiral EFT with pions and nucleons as the only explicit degrees of freedom. For this particular 
topology, the N'^LO corrections already provide the sub-subleading contribution, so that one can address the question 
of convergence of the chiral expansion. The pertinent results of our study can be summarized as follows. 



We worked out the N^'LO contributions to the 27r-exchange 3NF. The unitary ambiguity of the Hamilton 
operator can be parametrized at this order by three additional unitary transformations. We found that two 
of the corresponding "rotation angles", namely aio and an, are fixed in terms of the remaining one ag if one 
requires that the resulting 3NF matrix elements are finite (renormalizability constraint). The parameter ag 
does not enter the expressions for the 3NF at N^LO. These findings will impact the results for the remaining 
3NF contributions which arc not considered in this paper. 

In order to determine the low-energy constants q, di and e^ contributing to the 27r-exchange 3NF, we re- analyzed 
pion-nucleon scattering at order Q^ employing exactly the same power counting scheme as in the derivation of 
the nuclear forces. We used the available partial wave analyses of the pion-nucleon scattering data to determine 
all relevant LECs. The resulting values turn out to be rather stable and agree well with the determinations by 
other groups. 

With all LECs being fixed from pion-nucleon scattering as discussed above, we found a good/reasonable con- 
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FIG. 6: Chiral expansion of the functions A{q2) and B{q2) using the LECs determined from the fits to the KH iiN partial wave 
analysis. The dashed and solid lines are the same as in Fig. [S] The dotted and dashed-dotted lines show ^^^'(^2), B^^'{q2) 
with the LECs a taken from Eq. (|5.27p and (|4.25p . respectively. 

vergence of the chiral expansion for the functions A and B which parametrize the (static part of the) two-pion 
exchange 3NF. We also provide the recommended values of the LECs ci.34 which allow one to approximate the 
full, N'^LO result for the 27r-exchange 3NF by the N^LO expressions. 

As explained in the introduction, one generally expects important contributions to the nuclear forces associated with 
the intermediate A-excitations. The observed (reasonable) convergence pattern for the longest-range 3NF is not 
surprising given that effects of the A-isobar are, to a large extent, accounted for already in the leading contribution 
to this topology (i.e. at N^LO) through resonance saturation of the LECs C3.4. The situation is different for the 
two-pion-one-pion exchange and ring 3NF topologies, whose leading contributions at N'^LO are completely missing 
effects of the A-isobar. The corresponding N^LO corrections are, therefore, expected to be large and need to be 
worked out. Work along these lines is in progress. 
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Appendix A: Formal algebraic structure of the N*LO corrections 



In this Appendix we list the formal operator structure of the various N^LO contributions to the nuclear Hamiltonian 
relevant for the present calculations. A detailed discussion on the method of unitary transformation including the 
explicit form of the unitary operator at lower orders in the chiral expansion can be found in Rcf. [39| . 
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Here and in what follows, we adopt the notation of Refs. (la. [39l - l4l| . In particular, the subscripts a and h in H)^' refer 
to the number of the nucleon and pion fields, respectively, while the superscript k, gives the inverse mass dimension 
of the corresponding coupling constant^, see Ref. [30| for more details. The chiral order associated with a given 

contribution can easily be read off by adding together the dimensions k, of H^ ■ More precisely, it is given by 
^^ Ki — 2. In the above equations, 77 (A) denote projection operators onto the purely nucleonic (the remaining) part of 
the Fock space satisfying rf' = ij, X^ = X, rjX = Xij = and X + r] = 1. The superscript i of A* refers to the number of 
pious in the corresponding intermediate state. Further, E^^ denotes the total energy of the pious in the corresponding 



state, -Ett — X^i y^i^ +-^w> with li the corresponding pion momenta. Last but not least, the parameters ckg^io,ii 
parametrize the unitary ambiguity of the resulting nuclear Hamiltonian. 



Appendix B: Chiral expansion of the invariant tvN amplitudes g'^{uj,t) and h^{uj,t) 



In this Appendix we give the explicit expressions for the invariant amplitudes g {to, t) and h {to, t) which parametrize 
the pion-nucleon scattering matrix at first four orders in the chiral expansion. We use here the heavy-baryon approach 
with the nucleon mass being counted according to Q/m ~ Q^/A?^. This has the consequence that relativistic correc- 
tions are shifted to higher orders as compared to the standard approach based on the assignment m ~ A^. Apart 
from this difference, our results agree with the expressions given in Ref. [45.] (modulo one obvious misprint in that 
work) . 
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Contributions at order Q^: 
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Contributions at order Q^: 
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^ For 1/m-corrections, Ki corresponds to the inverse power of coupling constants plus twice the power of m ^. In particular, k = 2 for 
the nucleon kinetic energy term Hio- 



h+ 
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g\M^ {AglM^ {2M^ -t-2u;')+3u' {Mj - 2t)) g\ {^M^ - AM^t + t{t + 4a;^)) 

(Ji + J2) (4M2 uj^2tuj) 4d3 a;3 8 JgAf^ w disS^M^ (-2M2 + t + 2 u;^) 

Jo(l.) {g\ {Ml - C.2) (2Af2 - t -2uj^) + Qu:^) z^ oj^ - M^ {g\ {Mj ~ c.^) {2Ml ~ t -2u:^) + Qi^^) 

12F4a;2 + 967tF:^uj^ 

g\ (3Mg - 2 a;2) (2M2 - t - 2 ^2) _ ^2^ ^2 (13^2 + t) + u^^{t- 6Af 2) 

2887r2F4 w 

/20ft) ^ (-4 (2g^ + 1) Ml + Sgjt + t) {-AMJ + t + A^^) {g\ {2Ml - t + 2lj^) - 2uj^) 
12F* 16F2ma;2 

2{du-d^,)u; , 2Ji8g^M2 , .g^ JoH (c^^ - ^2) zg^ (^2 _ ^2)3/2 ^4(3^^2 + 4^2) ^^ (^ „ 4m2) 



F^ F^Lu QF^uj^ 487rF4w2 1447r2F4 

,3/2 



w 



^ zgi(a.2^Af2)^/- g^j^„ft)(4Af2_^) 4g^^Ml + 3glM^oj' gl{2Ml - t -2uj^) +2oj^ 

2471^4^2 + 8^4 96^i^4^2 + 8F2tow2 ' ^ ■'^> 

Contributions at order Q*: 
+ _ 2ci/2o(i)M2(M2-2<) / 72o(t) (4Af 4 _9Af 2^ + 2^2) 6M4 - 13M2t + 2f2 cj (-4Af2 + t + 4^2) \ 



C2 



F4 '" \^ 12F4 2887r2i^4 p2,^ 

C3/2o(0(2M4-5M2t + 2t2) Aeu{t-2Mlf 8ei5 l^^ (2Af g ^ t) IGeigc^^ 

2^4 + F2 + ^2 + ^2 

4M2 (2ei9 - 622 - 636 + 2r3Cif -^) (2M2 - t) 16(e2o + e35)Mgc^^ 8 (622 - 4e38 - Tadf -2) M4 

p2 + f2 + p2 

ici Ml uj y/uj'^-Ml ic2Uj^ y/uj^-Ml 
^ ^ 2^ 



g\Ml {2Ml - t - 2C.2) ^^c.2 - Mj {g\ {Ml - io^) (2Af2 - t - 2^2) + 6^4) 
'^^ 1^ 127rF4w 127rF4a; 

/^ ».g^ ( c^2 - m2) '/' (-2A/2 + t + 2 c^2) g\Ml {^2Ml + t + 2 u^) ^ tio \ 
64 



/i+ 



Ug\{u.^-Miy" g\Ml 

C3 - C4J - 



127rF4w 127rF4w 2i^^r 

3/2 



67ri^i w 67ri^4 , 



^_ ^ ^ / 4ffiJoH(M2-..2) -6(5ffi + l)M2 + 8gic.2 + t z.gj (Af 2 _ ^2) ^ ^2 _ ^2 
''''i^ 3i^4^ 1447r2F4 67ri^4^ 

J2o(t)(t-4Af2) ^\ e-17 (8Af 2 _ 4^) 8ei8C.2 8 (e2i - ^) Af2 

6F4 + F> j + F2 + F2 + F2 • ^ ■ ^ 

The low-energy constants and the kinematical variables entering these expressions are defined in sections |lTl IIIII and 
IIVI The loop functions are defined via: 






Kn it) = ■;= arctan — — - 

' Sny^ 2M^ 
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f m 1 y/l - 4M^/t , yi - 4M^/t + 1 
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